Liquid-Gas phase transition in Bose-Einstein Condensates 
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We study the effects of a repulsive three-body interaction on a system of trapped ultra-cold atoms 
in a Bose-Einstein condensed state. The corresponding s— wave non-linear Schrodinger equation is 
solved numerically and also by a variational approach. A first-order liquid-gas phase transition is 
observed for the condensed state up to a critical strength of the effective three-body force. 
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The experimental evidences of Bose-Einstein conden- 
sation (BEC) in magnetically trapped weakly interacting 
atoms 0-0] brought a considerable support to the theo- 
retical research on bosonic condensation. The nature of 
the effective atom-atom interaction determines the stabil- 
ity of the condensed state: the two-body pseudopoten- 
tial is repulsive for a positive s— wave atom-atom scat- 
tering length and it is attractive for a negative scattering 
length [Q. The ultra-cold trapped atoms with repulsive 
two-body interaction undergoes a Bose-Einstein phase- 
transition to a stable condensed state, in a number of 
cases found experimentally, as for 87 Rb Q, for 23 Na ||] 
and 7 Li However, a condensed state of atoms with 
negative s— wave atom-atom scattering_ length would be 
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unstable for a large number of atoms 

It was indeed observed in the 7 Li gas [3|, for which the 
s— wave scattering length is a — (—14.5 ±0.4) A, that the 
number of allowed atoms in the condensed state was lim- 
ited to a maximum value between 650 and 1300, which is 
consistent with the mean-field prediction ||. An earlier 
experiment suggested that the number of atoms in the 
condensate state was significantly larger than the theo- 
retical predictions with two-body pseudopotential. This 
is consistent with an addition of a repulsive three-body 
interaction, which can extend considerably the region of 
stability for the condensate even for a very weak three- 
body force. 

It was reported in Ref. || that a sufficiently dilute and 
cold bosonic gas exhibits similar three-body dynamics for 
both signs of the s— wave atom-atom scattering length 
and the long-range three-body interaction between neu- 
tral atoms is effectively repulsive for either sign of the 
scattering length. It was suggested that, for a large num- 
ber of bosons the three-body repulsion can overcome the 
two-body attraction, and a stable condensate will appear 
in the trap ||] . Singh and Rokhsar |l(| have also observed 
that above the critical value n (which is proportional to 
their — j c ) the only local minimum is a dense gas state, 
where the neglect of three-body collisions fails. 

In this work, using the mean-field approximation, we 
investigate the competition between the leading term of 
an attractive two-body interaction, which is originated 
from a negative two-atom s— wave scattering length, and 



a repulsive three-body interaction, which can happen in 
the Efimov limit jll| (\a\ — > oo) as discussed in Ref. || f\ 
We show that a kind of liquid-gas phase transition ap- 
pears inside the Bose condensate. 

The Ginzburg - Pitaevskii - Gross (GPG) nonlinear 



Schrodinger equation (NLSE) |16| is extended to include 
the effective potential coming from the three body inter- 
action and then solved numerically in the s— wave chan- 
nel. The dimensionless parameters are related to the two- 
body scattering length, the strength of the three-body 
interaction and the number of atoms in the condensed 
state. As particularly observed in Ref. |l7| ], to incorpo- 
rate all two-body scattering processes in such many par- 
ticle system, the two-body potential should be replaced 
by the many-body T— matrix. Usually, at very low en- 
ergies, this is approximated by the two-body scattering 
matrix, which is directly proportional to the scattering 
length a [|| . To obtain the desired equation, we first con- 
sider the effective Lagrangian density, which describes 
the condensed wave-function in the Hartree approxima- 
tion, implying in the GPG energy functional |T(| for the 
trial wave function "J: 
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In our description, the atomic trap is given by a rota- 
tionally symmetric harmonic potential, with angular fre- 
quency lo, and C\ gives the effective atom interactions up 
to three particles. 

The effective interaction Lagrangian density for ultra- 
low temperature bosonic atoms, including two [[| and 



*The physics of three-atoms in the Efimov limit is discussed 
in Ref. |12). This reference extends aprevious study of uni- 
versal aspects of the Efimov effect [g_3|. The relevance of 
three-body effects in BEC was also previously reported in 
Refs. (l^Jl^] , where it is discussed the stability of the numer- 
ical solutions. 
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three-body effective interaction at zero energy, is written 
as: 



m 61 



(2) 



where A3 is the strength of the three-body effective in- 
teraction and a the scattering length. 

The NLSE, which describes the condensed wave- 
function in the mean-field approximation, is variationally 
obtained from the effective Lagrangian given in Eq. (TP). 
By considering a stationary solution, ty(r,t) = e~ H ^ h 
where fi is the chemical potential and ip(r) is nor- 
malized to 1 and by rescaling the NLSE for the s— wave 
solution, we obtain 
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for a < 0, where x = y / zmujjh r and = 
iV 1 / 2 y / 87r[a[ rip(r). The dimensionless parameters, re- 
lated to the chemical potential and the three-body 
strength are, respectively, given by (3 = fi/Uuj and 173 = 
\3hujm 2 1 '{AttTi 2 "a) 2 . The normalization for <i>(a;) reads 
J dx\$(x) 1 2 — n where the reduced number n is related 
to the number of atoms N by n = 2N\a\yj2muj /%. The 
boundary conditions B in Eq.(^) are given by $(0) = 
and $0*0 — * Cexp(— a; 2 /4 + \J3 — |] ln(x)) when x — ► 00. 

The above equation, ([|) will be treated by numeri- 
cal procedures for non-linear differential equations, em- 
ploying the Runge-Kutta (RK) and shooting methods. 
However, it will be helpful first to consider a variational 
procedure fl8| , using a trial Gaussian wave-function for 
i/j(r) 
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where a is a dimensionless variational parameter. The 
corresponding root-mean-square radius is proportional to 
the variational parameter a, as ( r2 ) var = ct 2 3U/(2muj), 
while the central density is given by p c ,var(ct) — 
a~ 3 (mil) / Tifif^ 2 . The expression for the total variational 
energy is given by 



E var (a) = TiujN 
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In the same way, we can obtain the corresponding chem- 
ical potential , Eq. ( ||) : 
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The variational solutions of E var (a) are given, as a func- 
tion of n and #3 (where a < and 53 > 0), by finding the 
extrema of Eq. (||) with respect to a. 
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FIG. 1. In the lower part, we have a comparison between 
variational (solid curve) and exact (dashed curve) numerical 
calculations of the condensate energy as a function of the 
reduced number of atoms n for g-j, = 0.005 . In the upper 
frame we show five plots of the variational energy as a function 
of the variational parameter a for five particular values of 
n shown also in the lower frame. (I) (resp IV) corresponds 
to a small (large) n region where only one stable solution is 
encountered; (II) (resp III) to a small (large) n region where 
we observe three extrema for the energy; (C) corresponds to 
a particular n where we obtain two stable solutions with the 
same energy E\ = E2. E is given in units of (N%u))/n. 

In Fig. 1, we first illustrate the variational procedure 
considering an arbitrarily small three-body interaction, 
chosen as 33 = 0.005. In the upper part of the figure, 
we show five small plots for the total variational energy 
E, in terms of the variational width a. Each one of the 
small plots corresponds to particular values of n. For 
each number n we report the energy of the variational 
extrema in the lower part of figure 1. In region (I) where 
the number of atoms is still small, the attractive two 
body force dominates over the repulsive three-body force 
and just one minima of the energy as a function of the 
variational parameter a is found. That is also the case 
for $3 = 0. When the number of atoms is further in- 
creased (region (II)) two minima appear in the energy 
E (a) . An unstable maximum is also found between the 
two minima. The lower energy minimum is stable while 
the solution corresponding to the smaller a is metastable. 
This solution has a higher density and, consequently, its 
metastability is justified by the repulsive three-body force 
acting at higher densities. The minimum number n for 
the appearance of the metastable state is characterized 
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by an inflection point in the energy as a function of a. 
The value of n at the inflection point corresponds to the 
beak in the plot of extremum energy versus n because 
for larger n three variational solutions are found as de- 
picted in the lower part of figure 1. The attractive two- 
body and trap potentials dominate the condensed state 
in the low-density stable phase up to the crossing point 
(C). At this point, the denser metastable solution be- 
comes degenerate in energy with the lower-density stable 
solution and a first order phase transition takes place. 
Since the two solutions differ by their density this tran- 
sition is analogous to a gas-liquid phase transition for 
which the density difference between the liquid and the 
gas is the order parameter. In the variational calcula- 
tion this occurs at the transition number n ral.3 while 
the numerical solution of the NLSE gives 1.2. In region 
(III), we observe two local minima with different ener- 
gies, a higher-density stable point and a lower-density 
metastable point. The metastable solution disappears in 
the beak at the boundary between region (III) and (IV) . 
In regions (III) and (IV) the three-body repulsion sta- 
bilized a dense solution against the collapse induced by 
the two-body attraction. The qualitative features of the 
variational solution is clearly verified by the numerical 
solution of the NLSE, as shown by the dashed curve. 




FIG. 2. Central density p c , total energy E, chemical po- 
tential p, and average square radius (r 2 ^), as functions of the 
reduced number of atoms n. Several values of the dimension- 
less three body interaction strength g$ were used : gs = 
(short-dashed line), gz = 0.012 (dashed line), g$ — 0.015 
(dotted line), g 3 = 0.0183 (solid line), g 3 = 0.02 (dot-dashed 
line). The corresponding units are: (mu)/h)/(4:ir\a\) for p c , 
(Nhu>)/n for E, hcu for p, and %/(2muS) for (?" 2 ). 



In Fig. 2, considering several values of 173 (0, 0.012, 
0.015, 0.0183 and 0.02), using exact numerical calcula- 
tions, we present the evolution of some relevant physical 
quantities E, /i, p c and (r 2 ) as functions of the reduced 
number of atoms n. For (73 = 0, our calculation repro- 
duces the result presented in Ref. |q,|l7|l, with the maxi- 
mum number of atoms limited by n max « 1.62 (n is equal 
to \Clf\ of Ref. H). In the plot for the energy as a func- 
tion of n it is shown that for values of 53 > 0.0183 the 
phase transition is absent. At g 3 ps 0.0183 and n ps 1.8, 
the stable, metastable and unstable solutions come to be 
the same. This corresponds to a critical point associ- 
ated with a second order phase transition. At this point 
the derivatives of /i, p c and (r 2 ^ as a function of n all 
diverge. 

As shown in the figure, for < 53 < 0.0183, the 
density p c , the chemical potential p and the root-mean- 
squared radius (V 2 ) present back bendings typical of a 
first order phase transition. For each 173, the transition 
point given by the crossing point in the E versus n cor- 
responds to a Maxwell construction in the diagram of p, 
versus n. At this point an equilibrated condensate should 
undergo a phase transition from the branch extending to 
small n to the branch extending to large n. The system 
should never explore the back bending part of the dia- 
gram because as we have seen in figure 1 it is an unsta- 
ble extremum of the energy. From this figure it is clear 
that the first branch is associated with large radii, small 
densities and positive chemical potentials while the sec- 
ond branch presents a more compact configuration with 
a smaller radius a larger density and a negative chemical 
potential. This justify the term gas for the first one and 
liquid for the second one. However we want to stress that 
both solutions are quantum fluids. With 53 = 0.012 the 
gas phase happens for n < 1.64 and the liquid phase for 
n > 1.64. For g 3 > 0.0183 all the presented curves are 
well behaved and a single fluid phase is observed. We also 
checked that calculations with the variational expression 
of (r 2 ), p c and p are in good agreement with the ones de- 
picted in Fig. 2, following the same trend shown in Fig.l 
for the energy. 

Finally, in the lower frame of Fig. 3, we show the 
phase boundary separating the two phases in the plane 
defined by n and 33 and the critical point at n « 1.8 and 
<?3 ~ 0.0183. In the upper frame, we show the boundary 
of the forbidden region in the central density versus 53 
diagram. 

To summarize, our calculation presents, at the mean- 
field level, the consequences of a repulsive three-body ef- 
fective interaction for the Bose condensed wave-function, 
together with an attractive two-body interaction. A first- 
order liquid-gas phase-transition is observed for the con- 
densed state as soon as a small repulsive effective three- 
body force is introduced. In dimensionless units the crit- 
ical point is obtained when g$ w 0.0183 and n w 1.8. The 
characterization of the two-phases through their energies, 
chemical potentials, central densities and radius were also 
given for several values of the three-body parameter g 3 . 
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FIG. 3. Phase diagram of the Bose condensate. Central 
density p c in units of (muj/K) /(An\a\). 

The results presented in this paper can be relevant to 
determine a possible clear signature of the presence of re- 
pulsive three-body interactions in Bose condensed atoms. 
It points to a new type of phase transition between two 
Bose fluids. Because of the condensation of the atoms in 
a single wave-function this transition may present very 
peculiar fluctuations and correlations properties. As a 
consequence, it may fall into a different universality class 
than the standard liquid-gas phase transition, which are 
strongly affected by many-body correlations. This ques- 
tion certainly deserves further studies. 
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